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We calculate the effective resistivity of a macroscopically disordered two-dimensional conductor 
consisting of two components in a perpendicular magnetic field. When the two components have 
equal area fractions, we use a duality theorem to show that the magnetoresistance is nonsaturating 
and at high fields varies exactly linearly with magnetic field. At other compositions, an effective- 
medium calculation leads to a saturating magnetoresistance. We briefly discuss possible connections 
between these results and magnetoresistance measurements on heavily disordered chalcogenide semi- 
conductors. 
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The resistivity of most homogeneous materials (met- 
als or semiconductors) increases quadratically with mag- 
netic field H at low fields, and generally saturates at suf- 
ficiently large Exceptions may occur for materials 
with Fermi surfaces allowing open orbits, or for compen- 
sated homogeneous semiconductors, where the resistivity 
may increase without saturation, usually proportional to 
ff 2 i£. Under some special conditions, the magnetoresis- 
tance can be linear in magnetic field 3 . 

Recently, a remarkably large transverse magnetore- 
sistance (TMR) has been observed in the doped silver 
chalcogenides Ag2+aSe and Ag2+aTe 4 ^. In these mate- 
rials, over the temperature range from 4 to 300 K, the 
resistivity increases approximately linearly with H up to 
fields, applied perpendicular to the direction of current 
flow, as high as 60 T. Moreover, the TMR is especially 
large and most clearly linear at pressures where the Hall 
resistivity changes sigrA Because of this linearity, these 
materials may be useful as magnetic field sensors even at 
megagauss fields. 

But beyond the possible applications, the origin of the 
effect remains mysterious. According to conventional 
theories, such narrow gap semiconductors should have a 
saturating TMR. Furthermore, since these materials con- 
tain no magnetic moments, a spin-mediated mechanism 
seems unlikely. 

There are presently two proposed explanations for this 
quasilinear TMR. The first is a quantum theory of mag- 
netoresistance (MR)i. The second proposed mechanism-, 
is that this nonsaturating TMR arises from macroscopic 
sample inhomogeneities. Such inhomogeneities could 
produce large spatial fluctuations in the conductivity ten- 
sor and hence a large TMR, especially at large H. This 
explanation seems plausible because the chalcogenides 
probably have a granular microstructure^, and hence a 
spatially varying conductivity. 

The effective conductivity of media, with a spatially 
varying conductivity cr(x), has been studied since the 
time of Maxwell, but a relatively few studies have consid- 
ered the magnetoresistanco^ 10 ! 1 ^ 12 ' 13 ^ 4 ! 1 ^ 16 ! 17 ! 18 . For 
a three-dimensional medium, the TMR of an isotropic 
metal does indeed vary linearly in H, when a small vol- 



ume fraction p <C 1 of inclusions of a different carrier 
density is addecU^. But the TMR generally does not 
remain strictly linear at higher concentrations of p. If 
the inclusions are strictly insulating, then the TMR does 
remain asymptotically linear if the TMR is computed 
within the effective-medium approximation 1 ^, but its ex- 
act behavior is not known even in this case. Recent ex- 
periments on homogeneous semiconductors containing a 
gold inhomogeneity 20 show a hugely enhanced but not 
strictly linear room-temperature geometrical TMR (i. e. 
arising from inhomogeneities); this so-called extraordi- 
nary magnetoresistance has been successfully modeled, 
using finite-element techniques^!. 

The model of Refi& also assumes a film with a spa- 
tially varying conductivity. The inhomogeneities are de- 
scribed by an impedance network; the tensor nature of 
the magnetoconductivity is included by making each net- 
work element a four-terminal impedance. Their numer- 
ical solution suggests that, for the network to have a 
non-saturating TMR one needs (i) carriers of two dif- 
ferent signs, and (ii) a suitably defined average mobility 
(fx) r~j 0. When solved numerically and averaged over 
many disorder realizations, their model does indeed give 
a nonsaturating, approximately linear TMR over a broad 
field range. Obviously, it would be useful to have exact 
analytical statements to compare with these numerical 
results. 

In this Rapid Communication, we present an idealized 
model of a disordered semiconducting film in two dimen- 
sions. The model assumes a macroscopically inhomoge- 
neous film, consisting of two different types of conducting 
regions, denoted A and B, with areal fractions pa and 
Pb = 1 — Pa- In each region, the conductivity tensor is 
that of a Drude metal in a transverse magnetic field, but 
the density and the sign of charge carriers can be different 
in the two regions. We will show that, when pa = 1/2, 
and the charge carriers have opposite signs the TMR is 
asymptotically exactly linear at sufficiently strong mag- 
netic fields. Moreover, the linearity can extend down to 
quite low magnetic fields. The corresponding Hall coeffi- 
cient Ru,e is found to vanish, lip a ^ 1/2, the effective 
resistivity tensor p e cannot be calculated exactly. An ef- 
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fective medium approximation (EM A), which agrees with 
the exact result at pa = 1/2, predicts that the resistivity 
saturates for any pa ^ 1/2, and that Rn, e changes sign at 
PA = 1/2. All these results are in rough agreement with 
recent experiments^ [which are, however, carried out for 
three-dimensional (3D) samples; see below]. If the carri- 
ers have the same sign, no exact statements are possible, 
even at pa = 1/2. But even in this case the EMA pre- 
dicts a linear TMR precisely at pa = 1/2, though smaller 
than for carriers of opposite sign. 

We first prove the exact linearity of the TMR at pa = 
1/2 for carriers of opposite sign and opposite mobility, us- 
ing a duality argument. We consider a two-dimensional 
(2D) conductor with a spatially varying conductivity ten- 
sor <t(x), and denote the effective conductivity tensor by 
a e . <J e is a 2 x 2 tensor defined by (J) = cr e (E), where 
J and E are the position-dependent current density and 
electric field, and (...) denotes a spatial average in the 
limit of a large sample and suitable boundary conditions 
(as discussed, for example, in Refii^). a e is the quan- 
tity which would be measured as the sample conductiv- 
ity in an experiment. To calculate cr e , we use a duality 
theoremiS, which states that 



a e [a(x.)]a e [a 1 (x)} = I, 



(1) 



where I is the 2x2 unit matrix. Here cr e [cr(x)] denotes 
the effective conductivity tensor of a material whose local 
conductivity tensor is position-dependent and equal to 

(7(X). 

Thus, the product of a e for the system of interest, and 
that of a hypothetical "dual composite" whose local con- 
ductivity tensor Cd (x) is the local resistivity tensor of the 
original material, equals the unit tensor. 

We now apply this theorem to the following special 
case. Let the two components each have a free-electron 
conductivity, but carriers of opposite signs. For the first 
component 
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where <ja,o is the zero-field conductivity, the dimension- 
less magnetic field H — [iaB/c, where pa — cta/tra is 
an effective mobility of carriers of type A, niA is their 
effective mass, e > is the electron charge magnitude, 
and ta a characteristic relaxation time. For the second 
component, we assume 

: a B,yy = 7~u 0-2 ■ l 4 J 
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with the dimensionless constant k = — 1 (i. e., the two 
types of charge carriers have opposite signs). We also in- 
troduce [Lb — k[iA a s the effective mobility of type-B car- 
riers. Finally, we assume that the composite contains an 



areal fraction pi = 1/2 (i = A or B) of each component, 
and that the geometry is symmetric. " Symmetric" means 
that, if the components A and B were interchanged, o~ e 
of the film will remain the same in the thermodynamic 
limit. There are many geometries, both ordered (e. g. 
checkerboard) and random, which are symmetric by this 
definition. If we make the usual Drude assumption that 
0i,o = (i = A, B), where Uj is the number den- 

sity of carriers of type i, then eqs. 101-10) imply (i) that 
there are equal areal fractions of positive and negative 
charge carriers (but not that the total numbers of posi- 
tive and negative charge carriers are equal); and (ii) that 
the mobilities pa and ps are equal and opposite, so that 

(M) = Ei=A,B P^i = °; 

Given these assumptions, the tensors a a and o~b satisfy 
the remarkable relationship 



l + H 2 



(6) 



1/2 

where ctq = (ta.oCb.o) • Since we have an equal pro- 
portion of components A and B, distributed in some sym- 
metrical (and isotropic) fashion, the dual composite has 
a conductivity tensor 



ct(x), 



(7) 



where <x(x) means the conductivity of a composite in 
which the A and B components are interchanged. Since 
o~d is just a multiple of the original conductivity ten- 
sor <t(x), but with A and B components interchanged, 
and since by the assumption of a symmetric composite 
cr e [cr(x)] = <7 e [<j(x)], it follows that 



cr e [cr d (x)] = 



l + H 2 



Hx)]. 



(8) 



We now apply eq. (1) to this model, with the result 
l + H 2 



-o>(x)]=J. 



(9) 



A physically acceptable solution to cq. must have 
the diagonal elements of o~ e equal and positive, and off- 
diagonal elements equal and opposite. It is readily shown 
algebraically that the only such solution is 



cr e [cr(x)] = 



1 



- T a I. 



VT+lP 

The corresponding resistivity tensor p e is 



Pe 



Vl + H 2 I. 



(10) 



(11) 



The TMR is defined by the relation Ap SyXX (H) = 
[p e ,xx(H) ~ p e , xx(0)]/P e,xx(0)- For this model, 
Ap e ^ xx (H) = \fl + H 2 — 1 becomes linear in H for 
large enough H, and the corresponding Hall coefficient 
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FIG. 1: Calculated transverse resistivity p e , xx (H,PA) and 
Hall resistivity p e , xy (H,pA) (inset) for a two-dimensional 
model inhomogeneous semiconductor in a transverse mag- 
netic field, as calculated within the EMA for three different 
area fractions pa of component A. Both are given in units of 
1/cro = 1/ y/(TA,oo'B,o- The two components A and B have 
conductivities given by eqs. (2)-(5), with cta.o/cs.o = 2. The 
mobilities of the two carriers are assumed to have the same 
magnitudes: \(ia\ = |ms|- 

Rh = p X y{H)/H = 0. Thus, this calculation appears to 
reproduce the numerical results of Ref£, but analytically. 

Since the duality argument is not sufficient to deter- 
mine a e for pa ^ 1/2, we have used the EMA for such 
concentrations. The EMA is a simple mean-field approx- 
imation in which the local electric fields and currents are 
calculated as if a given region is surrounded by a suitably 
averaged environment. For the present model the EMA 
becomes^ 

PiScnil-TScTi)- 1 = 0. (12) 

i=A,B 

Here &<ii = <Ji — o~ e , and T is a suitable depolarization 
tensor. We assume that a a and <jb satisfy Oi^ xx = cri iVy ; 
&i,xy — —&%,yx- Then the components of er e satisfy 
<T e ,ra: = &e,yyi &e,xy = —^e,yx- We also assume that the 
two components A and B are distributed in compact, ap- 
proximately circular regions. Then Y = — I / (2a eiXX )2&. 
With these assumptions, eqs. (I12fl reduce to two coupled 
algebraic equations for 

o~ e,xx and o~ eX y which are easily 

solved numerically. 

To confirm that the EMA gives reasonable results, we 
have tested it for pa = Pb = \ , and a a and <jb given by 
eqs. with k = — 1. We find that the solution to 

the eq. (|12fl for the tensor a e is diagonal, and a multiple 
of the unit tensor; the diagonal elements are given by eq. 
pity. Thus, for pa = Pb, the EMA agrees with the exact 
duality arguments. 

To illustrate the EMA predictions for pa ^ \ , we cal- 
culate er e for &i given by eqs. The resulting el- 
ements of the resistivity tensor. p e ,xx — o~e,xx/[o~e xx + 

vlxy], Pe,xy = -^e,xy/[^l, xx + °l,xy\i are plotted ill Fig. 

1 for oao/cbq = 2. Evidently, and as can be shown ex- 
plicitly from the EMA equations, p exx is strictly linear 
in H only at pa — 1/2. At all other concentrations, 
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FIG. 2: Same as Fig. 1, but for pa = 1/2, and various choices 
of the mobility ratio k — psjpA- A positive or negative k 
means that the carriers have the same or opposite signs. 



Pe,xx{H) saturates (i. e. approaches a constant) at large 
H, but at a value much larger than p e ,xx(H — 0). It is 
easily shown that the saturation value of Ap e ^ X x(PA) = 
LiniH^oo [p e ,xx(H,p A )/pe(0,PA) - 1] oc 1/\pa - Pc\ on 
both sides of the percolation threshold p c = 1/2. Fig. 1 
also shows that the effective Hall resistivity p e . xy changes 
sign just at the concentration where p e ,xx varies asymp- 
totically linearly with H. 

We have also solved the EMA for a composite described 
by eqs. @ - © but for the more general case in which 
k 7^ — 1. Then k > and k < correspond respectively 
to carriers with mobilities of the same and opposite signs. 

In Fig. 2 we show the EMA results for this model. 
Specifically, we show p xx (H,pa) and p xy (H,pA) with 
Pa = 1/2, <ja,q = o'b.o, and several choices of k corre- 
sponding to carriers of both opposite and the same sign. 
The case k = 1 actually corresponds to a homogeneous 
free electron metal. For all other values of k, the TMR is 
asymptotically linear, the linear behavior is evident even 
at moderate fields (H ~ 1). However, the linear slope 
is larger when the carriers have opposite signs. We em- 
phasize that these results are obtained in the EMA. The 
duality arguments do not give any predictions for p xx ex- 
cept when the carriers have opposite signs and opposite 
mobilities. 

In Figs. 3 and 4, we plot the resistivity p xx and Hall 
coefficient Rh = Pxy/H as a function of pa for H = 1 
and H = 10. In both cases, we assume that ovlo = °\b,o 
and \pa\ = \pb\- Pxx has a peak at pa — 1/2, which 
sharpens, as a function of pa, as H increases. Similarly, 
the Hall coefficient Rh changes sign at pa — 1/2, and 
the change occurs over a narrower and narrower regime 
of pa as H increases. 

The present results agree qualitatively with the experi- 
ments of Lee et afe, which also show that the TMR peaks 
at pressures where the Hall coefficient changes sign. But 
this agreement should be viewed cautiously. In particu- 
lar, the measurements of Ref£ are carried out on a 3D 
sample, while our calculations are for a 2D system. The 
present work would also apply to a 3D system with a 
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FIG. 3: Transverse resistivity Pxx(H,Pa) and Hall coefficient 
Rh{H,pa) as a function of pa for H = 1, using the same 
model as in Fig. 1, with cjb.o = oa,o and pb = —pA- 
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FIG. 4: Same as Fig. 3, but for H = 10. 

columnar microstructure - that is, a system in which the 



conductivity tensor <t(x) is independent of the third di- 
mension, z - and the applied field B||z, but the samples 
of Ref.Si, if inhomogeneous, are most likely composed of 
small compact grains. We have calculated a e for a 3D 
granular sample with carriers of opposite signs, using 
the EMA, and find results similar to those shown here 
for 2D samples. These 3D calculations will be presented 
elsewhere^. 

The TMR of the present model is very large - 
Ap xx (H, 1/2) ~ 10 for H ~ 10 - and remains approx- 
imately linear down to fields as low as H ~ 1 — 2. By 
contrast, other models of TMR which arises from inho- 
mogeneities produce only a small TMR, or, if a large 
TMR, Ap xx (H) does not vary linearly with 75±£J£. 

In summary, we have presented a simple model of a 2D 
macroscopically inhomogeneous material, whose TMR is 
asymptotically linear in magnetic field, and whose corre- 
sponding Hall coefficient vanishes. The model has sev- 
eral unusual properties which make it likely to be real- 
ized only in special circumstances. First, eqs. ©-(HJ)) 
imply that the carriers have equal and opposite mobil- 
ities p. a = —p<B- Secondly, the linearity occurs only if 
the composite has a symmetric geometry at pa = 1/2. 
But given these features, the TMR, arising from a per- 
pendicular to the sample, is asymptotically exactly linear 
in B. To our knowledge, this is the only analytically solu- 
ble model for TMR due to macroscopic inhomogeneities, 
which produces a linear TMR at high concentrations of 
inhomogeneities. 
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